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Abstract 



Qh. The zero modes of the Dirac operator in the background of center vortex gauge 

field configurations in R 2 and R 4 are examined. If the net flux in D = 2 is larger 
than 1 we obtain normalizable zero modes which are mainly localized at the vortices. 



In D = 4 quasi-normalizable zero modes exist for intersecting flat vortex sheets 
with the Pontryagin index equal to 2. These zero modes are mainly localized at the 
vortex intersection points, which carry a topological charge of ±|. To circumvent 
the problem of normalizability the space-time manifold is chosen to be the (compact) 
torus T 2 and T 4 , respectively. According to the index theorem there are normalizable 
zero modes on T 2 if the net flux is non-zero. These zero modes are localized at the 
vortices. On T 4 zero modes exist for a non- vanishing Pontryagin index. As in R 
these zero modes are localized at the vortex intersection points. 
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1 Introduction 



There are two fundamental, principally non-perturbative phenomena of the strong inter- 
action which should be explained by QCD: confinement and chiral symmetry breaking. 
Soon after the advent of QCD chiral symmetry breaking has been explained by assuming 
that the QCD vacuum contains an ensemble of instantons and anti-instantons [[J. Lattice 
simulations have, however, shown that instantons cannot be responsible for confinement 
since they account for only about ten percent of the string tension 0, ^, |J. On the other 
hand, recent lattice calculations have produced mounting evidence that confinement is due 
to the condensation of center vortices in the Yang-Mills vacuum [5|, |6|. Since on the lattice 
the confinement phase transition is observed to occur at the same temperature at which 
chiral symmetry is restored, one expects that chiral symmetry breaking and confinement 
are triggered by the same mechanism. 

Since instantons cannot explain confinement one wonders whether center vortices are 
also capable of producing chiral symmetry breaking in addition to confinement. Indeed, 
it has been shown that in a Yang-Mills ensemble devoid of center vortices the relevant 
order parameter, the quark condensate vanishes |J. This does, however, not yet mean 
that center vortices produce chiral symmetry breaking in the QCD vacuum, since the 
field-configurations which produce it could be tied or attached to center vortices and si- 
multaneously removed with the latter. It is therefore still an open question whether center 
vortices produce chiral symmetry breaking and what the underlying mechanism behind 
this is. 

In order to get an idea how chiral symmetry breaking could be produced by center 
vortices let us first recall how it arises in the instanton picture of the QCD vacuum. In this 
picture the zero modes of quarks in the instanton background play a crucial role. These 
zero modes occur due to the topological charge of the instantons and are localized near 
their topological charge density |9], |10 |. 

In an instanton-anti-instanton ensemble the localized zero modes of the individual in- 
stantons spread out in space and form a continuum of quasi zero modes which by the 



Banks-Casher relation give rise to a non-zero quark condensate [|TT[] . In fact, independent 
of the instanton picture lattice calculations show a strong correlation between topological 
charge density and the quark condensate |T3| | . This correlation is mainly due to the zero 
modes which exist in topologically non-trivial gauge fields with non-vanishing Pontryagin 
index due to the Atiyah-Singer index theorem [THl and which are localized at topological 
charge. 

In center vortex field configurations topological charge is concentrated at the intersec- 
tion points and other singular points like twisting points Near these singular points 
we expect localization of quark zero modes which could play a similar role in the explana- 
tion of chiral symmetry breaking in the vortex picture as they do in the instanton picture. 
In a first attempt towards an understanding of chiral symmetry breaking in the vortex 
picture, in the present paper we study the quark zero modes in intersecting center vortex 
background fields. A more rigorous understanding of the mechanism of chiral symmetry 
breaking in the vortex picture is presently under investigation. 
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We examine zero modes of the Dirac operator in vortex backgrounds in two and four 
dimensions for different Euclidean space-time topologies. In sections [| and |3] fermionic 
zero modes in the background of vortices in M. 2 and intersecting flat vortices in M 4 are 
studied and their relation to solutions found by other authors [|16|, [T7|] is shortly discussed. 
In sections f|and [5] this analysis is repeated for space-time given by T 2 and T 4 , respectively. 



2 Fermionic zero modes in non-intersecting center vor- 
tex fields 

In D = 4 center vortices represent closed two-dimensional flux sheets. We are interested 
here in the quark modes in the background of such vortex sheets. Locally a center vortex 
sheet represents a two-dimensional plane. For simplicity we will consider in the following 
flat (planar) vortex sheets. In this case we have translational invariance parallel to the 
vortex sheets and the solution of the Dirac equation reduces to the two-dimensional problem 
in the plane defined by the 2 directions perpendicular to the vortex sheet. In this plane the 
vortex appears as intersection point. We will first study the Dirac equation in this plane, 
i.e. in D = 2. 

We will consider the zero modes of the Dirac operator in the background of Abelian 
gauge potentials representing Dirac strings and center vortices. We can consider these 
gauge potentials as living in the Cartan sub-algebra of a SU(2) gauge group. Having this 
in mind we say that a gauge potential describes a center vortex or Dirac string at a point 
Zq, if the Wilson loop around zq is —1 or +1, respectively. Equivalently, the magnetic flux 
carried by a center vortex or a Dirac string is given by $ 'center- — (m + 1/2) , m G Z or 
& Dirac = Tn , m G Z, respectivelyf]. 

The solution of the Dirac equation in two dimensions can be related to the theory of 
functions of a complex variable. We introduce the complex variable z by 

z = x + iy, z = x-iy, d z = - (d x - id y ) , d z = -(d x + id y ) (2.1) 
and a complex notation for the gauge potential. We define 1 

A z := - (A x - iAy) , A- z := ^ (A x + iA y ) = -~A~ Z , A x = 2i%A z , A y = 2i^A z , (2.2) 
where ^sA z is the imaginary part and 31A Z is the real part of A z . The Dirac equation 1 

nvJVA = AV, V=(^J, D lt = d lt + A li (2.3) 



in spinor components reads 



2i(d z + A z )ij 2 = A^i, (2.4) 
2i{d- z + A z )i h = A^ 2 . (2.5) 



^^Our conventions are summarized in Appendix 
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The most simple case is the free Dirac equation, i.e. A z = 0. Zero modes (A = 0) of the 
free Dirac equation are obviously given by analytic functions ip\ and anti-analytic functions 
ip2- A normalizable zero mode has to go to zero at infinity. But every (anti-) analytic 
function without singularities which goes to zero at infinity has to be zero everywhere. This 
means there are no smooth zero modes. In principle the zero mode may have singularities, 
in this case poles. But a pole in the zero mode is a non-integrable singularity and the zero 
mode would be not normalizable. Hence, there are no normalizable zero modes for the free 
Dirac equation. 

Next we consider the Dirac field in the background of a straight plane center vortex. 
Its gauge potential reads 



f(r 2 ) i _ f(r 2 ) i 

r 2 2 ' y X r 2 2 



A x = -y — — -, A y = x — — -, r = x +y = zz . (2.6) 



Here we introduced the profile function f(r 2 ) which fulfills /(0) = and f(r 2 ) — > 1 as 
r — > oo. The complexified gauge potential is then given by 

1 fir 2 ) f(r 2 ) 

A z := - (A x - iA y ) = ly. (-iy + x) = ?y- (z) . (2.7) 



0(x) : = / t^l&x' (2.1 



Introducing the function 

6(x) := I - 

2x 

the gauge potential A z can simply be written as 

A z = \d^(zz) . (2.9) 
Inserting the gauge potential into the Dirac equation we obtain the differential equations 

2i(d s + ±d s <j>(zz))ik = A^i, (2.10) 

2i(d- z -^d- z (j>{zz)\^x = A^ 2 . (2.11) 

We are mainly interested in the zero modes A = 0. Furthermore, let us first consider 
an idealized vortex with f(r 2 ) = 1 (this function obviously does not have the properties 
of a profile function, because /(0) = 1 ^ 0). For f{r 2 ) = 1 the function becomes 
4>{r 2 ) = |log(r 2 ), which is - up to a factor - the Greens function of the Laplace operator 
in two dimensions. For A = and 4>(r 2 ) = | log(r 2 ) the differential equations (|2.10| , |2.11|) 
can again be easily solved 

= {<fzl) X i{z), (2.12) 
^ 2 = (Vz^^Mzj, (2.13) 



where xi an d Xi are analytic functions of z and X2(z) is the complex conjugate of Xz( z )- 
Choosing e.g. xi = 1 we get — "sA 7 ; i- e - a real-valued function of x and y. The zero 
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modes (|2.12| , |2.13|) of the Dirac operator in the background of a single center vortex are 



not normalizable. This is because normalizable analytic functions X1/2 have to approach 
zero at infinity, and therefore, have to be identically zero or have a pole which yields a 
non-integrable singularity as was already discussed above. 

This result is in accord with the index theorem for U(l) gauge fields in D = 2. The 
theorem states that the difference between the numbers of right- and left-handed fermionic 
zero modes, say n, in the background of a U(l) gauge potential on M 2 is related to the 
total flux $ = ^ J R2 F by @ 

n = [$] , (2.14) 

where [x] is the largest integer smaller than x G K. (i.e. [1] = 0). In the case of a single 
center vortex on M. 2 the total flux is equal to 1/2, i.e. the number of left-handed zero 
modes is equal to the number of right-handed zero modes (in the present case there are no 
normalizable left- or right-handed zero modes). There exists an extensive literature about 
the Dirac operator in the background of Abelian gauge fields in two dimensions. A few 
remarks about some papers concerning also eigenmodes of such Dirac operators are in order 



here. In [16j the authors considered zero modes of the Dirac operator in the background 
of a flux tube on a disc with finite radius and APS boundary conditions. In this case the 
index n of the Dirac operator receives contributions from the boundary of the space-time 
manifold]]. In [17| self-adjoint extensions of the Dirac operator in the background of a 



singular magnetic vortex have been examined. Due to the boundary conditions used in 
these papers the encountered eigenmodes differ from the ones found here. In our opinion it 
is not clear which boundary conditions are the correct ones in a physical context. Therefore 
we will consider a smooth gauge potential representing a smeared out vortex. This way 
we get rid of the singularity of the vortex field and avoid the discussion of self-adjoint 
extensions of the Dirac operator at the position of the vortex. We choose a profile function 
/(r 2 ) = r 2 /(e 2 + r 2 ) , eel, i.e. we work with the gauge potential 



r 2 + e 2 2 r 2 + e 2 2 



For this profile function / we obtain <fr(r 2 ) = | log(r 2 + e 2 ). For A = the differential 
equations (|2.10| , |2.11 ) can be solved also in this case: 



fa = + Xi(z), (2.16) 

V> 2 = fv/HTi 2 )' 1 ^). (2.17) 



For X1/2 — 1 these solutions are obviously not normalizable. With the same arguments as 
before we conclude that there are no normalizable zero modes. 

As discussed already above in connection with the index theorem the asymptotic be- 
havior of the zero mode of the Dirac operator changes if we change the magnitude of the 



2 We consider in the present paper space-time manifolds with the topology M , S N or T N . The latter 
two have no boundary whereas on $L N the gauge potential vanishes at infinity. Therefore in our case there 
are no boundary contributions to the index. 
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magnetic flux. Multiplying the gauge potential A in ( |2.15| ) by a factor p e R + the flux $ 
becomes p/2 and the solutions VV 2 > c ^ e< l s - ( liH ' |2 • 1 7| ) , change to (ipi/2) p - This means 



that for p > 2 we get a total flux $ > 1 and a normalizable zero mode with chirality —1: 



^ 2 = (^+£2)-', Vi = 0. (2.18) 

The probability density of this zero mode has a maximum at r = 0, i.e. at the location of 
the vortex. 

Since the differential equations ( p.lOj , |2.11| ) are of first order type it is also possible to 



explicitly write down zero modes for multi-vortex background fields, because the multi- 
vortex gauge potential can be written as the sum over the gauge potentials of several 
simple vortices. The solution of the corresponding multi-vortex differential equation is 
then simply the product of the solutions to the several one-vortex differential equations. 
As an example take two vortices, one at z = a and a second at z = b. The solutions ipx/2 
then read 

Vi = -a)(z-a)+e*tf(z- b)(z - b) + e^j X i(z) , (2.19) 

^2 = (V(z -a)(z-a)+e*tf(z- b)(z - 6) + e^j X&j ■ (2.20) 

There is one quasi-normalizablePI zero mode: (O,^)- This zero mode is obviously 
localized at the centers of the vortices (at z = a and z = b) and on the line between them, 
cf. fig. |l[ A similar calculation on the compact sphere S 2 yields a normalizable zero mode, 
cf. Appendix fB|. According to the index theorem we obtain normalizable zero modes as 
soon as the net flux of the vortices exceeds 1 (i.e. two units of the flux of a center vortex). 

3 Fermionic zero modes of the Dirac operator for in- 
tersecting center vortex fields 

So far we have considered the fermionic modes in the background of parallel non-intersecting 
flat center vortex sheets. Due to the translational invariance parallel to the vortex sheets 
it was sufficient to study the Dirac operator in the plane perpendicular to the vortex sheet, 
where the vortices appear as intersection points. In the confined phase the center vortices 
percolate ||19|| . The vortices then have arbitrary directions and also intersect. In intersec- 



tion points all four space directions participate and, consequently to study the fermionic 
modes in the background of intersecting vortices we have to consider the full 4-dimensional 
Dirac operator. 

In the following we consider four vortex sheets consisting of two orthogonal pairs of 
parallel vortex sheets. One vortex pair is given by two planes parallel to the x\ — X2 plane 
located at = and X3 = ±0. The other vortex pair consists of two planes parallel to 



'Strictly speaking we have to increase the magnetic flux infinitesimally to get a normalizable zero mode. 
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Figure 1: Probability density of the 
zero mode in the background of two 
center vortices in D = 2 (see eqs. ( \2.l9i , 
TWj )) for e = 0.01, a = 1 and b= -1. 



Figure 2: Probability density of the 
zero mode in the background of four 
intersecting vortex sheets is shown in 
the subspace £4 = 22 = for a = b = 1 
and e = 0.01. 



the £3 — £4 plane located at £2 = and x\ = ±b. The four vortices intersect in four points 
given by £4 = £2 = 0, £3 = ±a, X\ = ±b. Each of the intersection points carries topological 
charge ±1/2. If we choose the orientation of the flux of the two vortex sheets to be equal 
we find a total Pontryagin index of v = ±2, otherwise (if the fluxes are anti-parallelly 
oriented) the Pontryagin index vanishes. For definiteness we choose the direction of the 
flux in parallel vortex pairs to be the equal. Then the gauge potential can be chosen as 



At 

A 3 
where 



-£ 2 - 



%2- 



—£ 4 - 



£4- 



2 \i 



1 

2' 



.4, 



Aa 



(£1 + b) 



(£3 + a. 



fist 



+ (£1 - b) 



+ (£ 3 - a 



(£ 3 ± a) 2 + x 
r 2 /(r 2 + e 2 ). 



2 

■1 1 



[x\ ± b) 2 + x\ , a, 6 G 



The field strength F^ u of the above gauge potential is given by 

2s 2 2e 2 



F\2 



[sl + e 2 ] 2 
2e 2 



[ri + e 



2\2 



+ 



+ 



[sl + e 2 ) 2 
2e 2 



(r 2 _ + e 



2\2 



Splitting the field strength into (anti-)self-dual components 



f(s 2 -)} 


i 


si I 


2' 


f(r 2 _y 


i 




2' 



(3.21) 
(3.22) 



(3.23) 
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where 77" and 77" are the t' Hooft symbols and a = F\ 2 + F 34 , (3 = iq 2 — -F34, one finds 
that in the limit e — > the field strength is self-dual (/3 = 0) at the intersection points. 

To determine the zero modes of the Dirac operator in the background of the considered 
intersecting vortices we introduce complex variables u and v [EUf 



u = £3 + ix^ , v = X\ + ix 2 
and the corresponding complex derivatives d u and d v 



1 (d X3 - id X4 ) , d v = - (d xi - id X2 



(3.24) 



(3.25) 



We also introduce two complexified gauge potentials A u and A v as in equation (|2.7| ) 
These gauge potentials can be written as 



where 



.4, 



log 



.4, 



(3.26) 



((u — a)(u — a) + £ 2 )((u + a) (u + a) + e 2 ) 
{{v -b)(u-b) + e 2 )((u + b)(u + b)+ e 2 ) 

The Dirac operator in the background of the gauge potential (|3.21| ) reads 

/ iD 4 + D 3 Di - iD 2 \ ( Vi \ 

n n 



(3.27) 





D!+iD 2 iD 4 -D 3 

iD 4 + D 3 Di - iD 2 

\ Di + iD 2 -iD 4 -D 3 



^2 
^3 



(3.28) 



where D M = <9 M + v4 M is the covariant derivative. For the considered vortex gauge poten- 
tial ( p.21| ) the eigenvalue equation ( |2.3j ) for the Dirac operator reduces to four coupled 
differential equations: 









A^2 






A^3 




{ 2d u 


A^4 




( 2d- v 



du<f> ) fa + i[ 2d v + -d v (j> ) ip 4 , 



dv<p ) ^3 + M -2<9 U - -<9 U ) ^ 4 , 



1 

2' 



which can be easily solved for A = 0, since in this case the upper and lower components of 
the Dirac spinor decouple. The A = solutions read 



^( r 2 +£ 2 ){r 2_ +£ 2 ) \ y ( ^ +£2)(s 2 +£ 2) Xl ^ v) 



= ^j(r\ + e 2 )(r 2 _ + e 2 ) (^j{s\ + e 2 )(s 2 _ + e 2 )) X 2(u,v) 



^3 = i/(r 2 + + e 2 ){r 2 _ + e 2 ) $j ' {s\ + 5 2 )(s 2 _ + e 2 ) X3 (u, v) , 

0>4 = (d{r\ + e 2 ){rl + e 2 )) (t/tf. + e 2 )(s 2 _ + e 2 )^ Xi (u,v) . 

The analytic functions Xi have to be chosen constant to avoid non-integrable singularities 
as was discussed for the D = 2 case. The only spinor component going to zero at infinity is 
■04. Therefore, only the spinor with non- vanishing component 04 (and all other components 
zero) yields a normalizable zero mode (normalizable up to a logarithmic divergence from 
the integration r± —>■ oo and s± — > oo - but these divergences become integrable if the 
vortex gauge potential is multiplied with p = l+ a,a>0). Embedding the U(l) gauge 
group into an SU(2) gauge group results in a second zero mode with opposite iso-spin 
(corresponding to inverting the sign in front of the vortex gauge potential (|3.21 )) and the 
same chirality as the above zero mode. This second zero mode is given byQ 

01 = 02 = 04 = , 

-1 



03 = ^ (r 2 + + E 2 )(r 2 _ + E 2 )j f ^(4+ £ 2) (s 2 +£ 2) 

The probability density of the zero modes is peaked at the four vortex intersection 
points and at the vortex sheets, cf. fig. H 



4 Center vortices and Dirac equation on the 2-torus 

In this chapter we consider zero modes of the Dirac operator in the background of Abelian 
gauge potentials representing Dirac strings and center vortices on the torus T 2 . There is 
a variety of reasons for studying T 2 in addition to IR 2 : first, T 2 allows to use the Atiyah- 
Singer index theorem in a stringent fashion and allows for normalizable spinors in the 
background of integer flux. Second, from a physical point of view one would want the 
mechanism of chiral symmetry breaking to depend on local quantities rather than on global 
characteristics like boundary conditions imposed from the manifold. However, individual 
solutions of differential equations may depend strongly on the boundary conditions. So if 
the envisaged mechanism of chiral symmetry breaking contains a truncation to a subset 
of modes of the Dirac operator, it is not a priori clear that the mechanism is indeed 
independent of the boundary conditions. Therefore the proposed mechanism has to be 
checked for different topologies of the space-time manifold. Third, the torus simulates 
a periodic arrangement of vortices. This is much closer to a percolated vortex cluster 
than a single vortex in R 2 . It will turn out that the zero modes are again localized at 
the position of the vortex, thus strengthening the point of view that the zero modes are 

4 Embedding the U(l) gauge group into SU(2) gives the Pontryagin index 2 for this configuration. This 
corresponds to the 2 observed zero modes. 
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indeed influenced by local properties of the gauge potential only. Finally, the torus is the 
space-time manifold that is also used in lattice calculations. 

At this point it should be mentioned that zero modes of the Dirac operator in the 
background of ideal vortices on the torus are also relevant to find instantons on a (dual) 



torus by the Nahm transformation [23, 24 



4.1 Periodicity properties of the gauge potential on the torus 

We first consider flat vortices so that we can restrict ourselves to a two-torus where the 
vortices appear as piercing points. Instead of working on the compact torus we can work on 
its covering manifold - on T 2 = R 2 = C. But in this case we have to demand that physical, 
i.e. gauge invariant, observables are periodic. This implies that the gauge potential has to 
be periodic up to gauge transformations. We choose the two-torus to have circumferences 
1 and t. Then the gauge potential A has to fulfill "quasi-periodicity" conditions 

A{z + 1) = A u * (z \z), A{z + it) = A u ^ z \z) , (4.29) 

where A v = f/ _1 AJJ + U~ 1 d tl U is the gauge transform of A. The transition functions 



U x , U y have to fulfill the cocycle condition pi 



U x (z)U y {z + 1) = U y {z)U x {z + it) . (4.30) 

We can always choose U x = 1 and U y independent of y. This implies periodicity of A(z) 
in x-direction. Furthermore, U y as function of x defines a mapping from S l into U(l). 
Such mappings fall into homotopy classes tti(U(1)) = n^S 1 ) which are characterized by a 
winding number n. A simple calculation shows that this winding number is related to the 
magnetic flux $ through the torus 

$ = -n, (4.31) 
reflecting the quantization of the magnetic flux through the torus. 



4.2 Dirac string on the torus 

Before we write down the gauge potential of center vortices on the torus let us consider 
the gauge potential of a single Dirac string. The reason is that a single center vortex does 
not exist on a torus while a single Dirac string does. 

On the two-torus the gauge potential of a singular point-like object can be expressed 
by means of the theta functions 

6(z, it) = e-" T " 2+2 ™ 2 , teR + (4.32) 

which are analytic in z and obey the periodicity properties 

9(z + l,iT)=9(z,iT), 6(z + iT,iT) = e 1TT - 2mz 6(z,tT). (4.33) 
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The only zeros of this function are at the points f22] 

z=(m+l/2) + (n + l/2)ir, m,neZ. (4.34) 
For subsequent consideration let us also introduce the real-valued function [^3|, Q 

(/) (z, zq) = - log (e T (z, z )9 T (z, zofj , 9 T (z, z ) := 9{z + 1/2 + l/2ir - z , vr) , (4.35) 

which (up to a factor of 27r) represents a (non-periodic) Greens function of the Laplacian 
on the two-torusQ. By Taylor expanding 9 T (z, zq) around its zero at z = z one finds that 
the Greens function <f>o(z, zq) behaves near zq as log|z — zq\. By means of this Greens 
function the gauge potential of a Dirac string on the two-torus can be expressed as 



A, 



-id,, 



A q . — id x 



Using the same notation as in equation Q2.2|) 

A z : = 1/2(4, - iAy) , A x = 2iQ(A 

we obtain 

I 

-id y (po + d x (f) ) = d z (f) 



A ? 



1 

2' 



A y = 2i$l{A z 



ld z 6 T (z,z ) 



(4.36) 



(4.37) 



(4.38) 



2 e T {z,z ) 

where we used that 9 T (z, zq) is an analytic function of z and 6 T (z,zo) is an anti-analytic 
function, i.e. it does not depend on z. The periodicity properties of A Z ,A X and A y can 
simply be derived using equation (|4.33| ) 



A z (z + 1) 
A y (z + 1) 
A x (z + 1) 



A z (z), A z (z + ir) = -iri + A z (z) , 

A y (z + ir) = A y (z) , 

A x (z), A x (z + ir) = -2m + A x (z) . 



(4.39) 
(4.40) 
(4.41) 



This means that the gauge potential fulfills equation (|4.29| ) with U x (x, y) = t and U y (x, y) = 
exp (—2irix)l. Furthermore, computing explicitly the flux going through the torus, we ob- 
tain (using Stokes theorem) 



(i,o) 

— I / A x dx + 



1 

2~i 



1 

2~TTi 



(0,0) 

1 



(l,r) 



A y dy + 



(1.0) 



(0,r) 



A x dx + 



(1,t) 



(0,0) 



Aydy 



(0,r) 



(A x (x,0)-A x (x,r))dx+ / (A y (l,y)-Ay(0,y))dy =1, (4.42) 



which is consistent with ( |4.31|) . The Dirac string is located at the point where 9 T (z, Zq) = 0, 
i.e. at z = Zq + m + inr , m,n G Z. The field strength of the Dirac string configuration 
vanishes at points where 9 T (z,z ) ^ 0, because 

F xy = d x A y - d y A x = i(d 2 x + d 2 y )<j) = R(4d g d g <i>) = 2W- z {d z 9 T {^ z )/9 T (z, z )) = . (4.43) 



'On the torus there is no periodic Greens function of the Laplacian 
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Here we used again that 9(z, zq) is an analytic function, i.e. it is independent of z. On 



the other hand the flux through the torus is 1, cf. ([4.42 ), from which we conclude that we 
have a Dirac string (which is represented by a point in D = 2) at the zero of the function 
9 T (z, Zq), i.e. at z = z + m + nir. 

The Dirac string can also be written as a pure gauge. If we define the U(l) gauge 
function 

which is singular at the zeros of 9 T (z, z ), the gauge potential 

A, = g-%9 (4.45) 

becomes 



A z = l(A 9 -iA v ) = g-%g = jy^dJtkdL (4.46) 



1 d z 9 T (z,z Q ) 



(4.47) 



2 T {z,z o ) 

which agrees with equation ( [4.381) . The periodicity properties of g are as follows: 

g(z + 1) = g(z) , g(z + it) = -e~ w g ( z ) . (4.48) 

Since the gauge potential of the Dirac string is a pure gauge we can simply write down the 
zero modes of the corresponding Dirac operator. The zero mode is nothing but the gauge 
transformation of a constant Dirac field. 

In the following we will consider only smeared out vortices on the torus. This way we 
avoid the discussion of boundary conditions of the the spinor field at the position of the 
singularity of an ideal vortex. This discussion would involve the problem of self adjoint 



extension of the naive Dirac operator, see e.g. |16, [17], [25 . 

The gauge potential of a smeared out Dirac string at the point Zq is given by equation 
( |4.36| ) with the function <j>o(z, Zq) replaced by 



<j>(z,Zo) = -\og(9t(z,z )9t(z,z )+6;(z,z )6-(z,z )), (4.49) 
9${z,z ) = 9(z + ^ + ^T-z ±e,ir), eeR+, (4.50) 



i.e. by 



A z = l/2(A x -iA y )=d z <f>(z,z k ) (4.51) 

(4.52) 



1 9+(z, z k )d z 9+(z, z k ) + e~(z, z k )d z 9 T (z, z k ) 



2 9+(z,z k )9+(z,z k ) + 9-(z,z k )9-(z,z k ) 
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By Taylor expanding around z k we obtain the behavior of the gauge potential for small 
distances r from the center zq of the Dirac string: 



A, 



-V-. 



A, 



x- 



(r 2 + e 2 ) ' y (r 2 + e 2 ) 

This gauge potential indeed represents a smeared out Dirac string. In the limit e 
becomes the gauge potential of an ideal Dirac string |Ll| on the torus. 



(4.53) 



4.3 Fermionic zero modes of center vortices on the torus 

A center vortex living in the Cartan sub-algebra can be represented by half the gauge 
potential of a Dirac string. However, for a single center vortex it is not possible to relate 
the gauge potential at z with the gauge potential at z + ir by a transition function U y 
which is periodic in x. Instead we need an even number of center vortices on the torus, in 
accord with the quantization of magnetic flux through the torus, cf. ( 4.31 ). We consider 
the configuration, see eqs. (|4.49| , f4.50|) , 

A z = ^(d z cf ) (z,z 1 )+d z <j ) (z : z 2 )) (4.54) 



1 9+(z, z k )d z 6+(z, z k ) + 6 T (z, z k )d z 6 T (z, z k ) 
A hi z k )6+(z, z k ) + 6-(z, z k )6-(z, z k ) ' 

which consists of two (smeared out) center vortices at the points z\ and Z2 and fulfills the 
periodicity properties (|4.29|) with transition functions U x — 1 and U y = exp (— 2irix)t = 
exp (— ni{z + z))t. 

We are interested in the A = eigenfunctions of the Dirac operator in the background 
of these two vortices. In accord with the boundary conditions to the gauge field ( |4.29 ) we 
impose the boundary conditions 

ip(z + 1) = U^z)- 1 ^) = tjj(z) , ip(z + ir) = U y (z)- l i)(z) = e ni(s+z) ^(z) (4.55) 

to the Dirac spinor if>. The Dirac eigenvalue equation for the two spinor components reads 

i(2d z + 2A z )^ 2 = Xip u (4.56) 
1(2^ + 24^1 = A^ 2 , (4.57) 

where A z = 1/2(A X + iA y ) = —A z . Because of the simple form of the two- vortex gauge 
potential (|4.54p the zero modes of the Dirac operator can be found explicitly. A short 
calculation yields 

/ 2 \ 1/4 

^ = ( II (et(z,Zk)enz,zk) + e;{z,z k )e;(z,z k )U (4.58) 
^2 = l]l(et{z,z k )e+{z,z k ) + e-{z,z k )e-(z,z k ))\ x 2 (^). (4.59) 
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The analytic functions Xi( z ) an d X2{z) have to be chosen such that ipi/2 are normalizable 
and fulfill the periodicity properties (|4.55|) . To render ip\/2 normalizable the functions X1/2 
must not have poles^. Inserting equations (|4.59 , 4.58 ) into ( 4.55Q and using the periodicity 



properties ( |4.33| ) of the theta function we arrive at 

X x(z + ir) = e (-^ +2 < z+ ^-^ + ^)) X i(z), Xi(z+ 1) = , (4-60) 
X 2(z + ir) = e ^ 2 < z+ ^-^ +z ^)x2(z), + 1) = ;&(*) • (4.61) 

If we require analyticity of X2 on the whole torus then this function is fixed (up to a factor) 
by the periodicity properties ( [4.61| ) to be the theta function 



X2 (z) = 9(z + l -r- + z 2 ), it) . (4.62) 



This is proven in Appendix 0. The function ip2, cf. ( [4.59| ), is obviously normalizable, be- 
cause it has no singularities. The only zeros of X2{z) are at the points z = ^+^$s(zi + z 2 ) + 
m + nir. On the other hand the required periodicity properties for xi ( |4-60| ) show that 
there is no (non-trivial) function Xi which is analytic on the whole torus. The function xi 
has to have at least one pole. Such a solution is given by the inverse of a theta function 

Xx{z) = + l -r - % -% Zl + z 2 ),ir) . (4.63) 

The poles of this function are at the points z = z mn — | + i£s{z\ + z 2 ) +m + nir , m,n G Z. 
Our considerations show that there is only one normalizable zero mode (with %pi = 0). If 
the component ipi of the Dirac spinor is not identically zero then it would have a pole at 
z = z mn , because the pre-factor of xi m equation (|4.58| ) is non-zero on the hole torus. This 
pole would yield a logarithmic divergence. 

In Appendix multi-vortex configurations and the corresponding zero modes of the 
Dirac operator are presented. 

There is another interesting point related to the periodicity properties given by ([4.29 



4.55| ). Multiplying the transition functions U x and U y by constant phases, say e la and e 1/3 , 



resp., the center vortex and Dirac string gauge potentials, (|4.54|) and fl4.36[) , still fulfill the 
periodicity properties with the new transition functions U x = e~ ia t and U y = e - 2mx - t i 3 f_ 
But these new periodicity properties change the zero modes of the Dirac operator by 
changing the periodicity properties of the analytic functions Xi/2(z). The new solution 
reads 

Uz) = e—6(z + i r-^r- Igfa + z 2 ) + JL, ir) (4.64) 

I 111 I Z7T 

The zeros of this function are at the points z = h + i$$(zi + z 2 ) + ^-t— f- +m + nir, i.e. the 
zeros are shifted by %OLT ^ compared to the original case, where U x = 1 and U y = e~ 2mx l. 



A pole of X1/2 would yield a non-integrable (logarithmically divergent) singularity in the norm of the 



spinor ip, because the pre-factors of X1/2 m equations (4.58,4.59) are non-zero on the whole torus, since 
we only consider smeared out vortices. 
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Figure 3: Probability density of the 
zero mode in the background of two 
center vortices on T 2 for r = 1, e = 
0.01, z l = 0.25i and z 2 = 0.751 



Figure 4: Probability density of the 
zero mode in the background of four in- 
tersecting vortex sheets on T 4 shown in 
the subspace x\ = x 3 = for u\ — V\ = 
0.25z , u 2 = v 2 = 0.75z and e = 0.01. 



The change of the transition functions by multiplying with constant phases e~ ia and 
e -i/3 j g equivalent to introducing a constant background gauge potential A z = |(~ + ia) 
and leaving the transition functions unchanged. 

The above considerations have shown that for a two- vortex gauge potential (smeared out 
vortices) we get exactly one normalizable zero mode which has exactly one zero on the torus. 
The position of the zero depends on the imaginary parts (y-coordinates) of the positions 
of the vortices and on the periodicity properties of the gauge potential (or equivalently 
on the presence of a constant background field) and of the spinor field. Furthermore, the 
probability density of the spinor field is peaked at the positions of the vortices. 

5 Fermionic zero modes for intersecting vortices on 
the 4-torus 

As in the case of space-time M 4 the fermionic zero modes for intersecting vortices on T 4 
can be explicitly written down. 

We consider four smeared out center vortex sheets consisting of two orthogonal pairs of 
parallel vortex sheets intersecting in 4 points as in section ^. Introducing complex variables 
u and v, cf. eq. ( |3.25| ), the complexified gauge potential A z can be chosen as 

A z = d u (<f)(u,u 1 )+<f)(u,U2)) + d v (<j)(v,vi)+(f)(v,V2)), 



where the function (f>(z, Zk) is defined in equation ( |4.49|) . As in the case of the space-time 



manifold M 4 there is only one normalizable zero mode given by 

Ipl = 1p2 = ^3 = , 
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ip4 = y\\\6^(u,u k )6+(u,u k ) + 6 T (u,u k )6 T (u,u k ) 

Y[ (o?(v,v k )d+(v,v k )+0~(v,v k )0-(v,v h ))\ Xi(u,v) 

vfc=l 

where Xi{ u i v ) is an analytic function of u and v given by 

Xa(u, v) = 0(u + -r - 2^( M i + ^2), ir)0(v + -r - -O^ + v 2 ), it) . (5.65) 

The probability distribution of this zero mode in the plane 9fai = §lv = for u% = u 2 = 
0.25i , m 2 = v 2 = 0.75i is plotted in fig. [|. 

6 Concluding remarks 

In the present paper we have studied the properties of fermionic zero modes in a center 
vortex background field. We have demonstrated that these zero modes are concentrated 
at the localization of the center vortices. In accord with this the probability density of 
these zero modes is sharply peaked at the vortex intersection points which carry (localized) 
topological charge 1/2. This result is consistent with the localization of the fermionic zero 
modes in an instanton background field at the instanton center. In fact lattice calculations 
show a strong correlation between the topological charge density distribution and the 
distribution of the quark condensate {q(x)q(x)) , the order parameter of chiral symmetry 
breaking. Given the localization of the quark zero modes at the localization of topological 
charge, we expect the quark zero modes in the vortex background field to play a crucial 
role for the spontaneous breaking of chiral symmetry in the vortex picture. This will be 
subject to future investigations. 
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A Conventions 



We choose the generators of the gauge group to be anti-hermitian. Therefore the com- 
ponents of the gauge potential are anti-hermitian, e.g. purely imaginary for the gauge 
group 2/(1). The magnetic flux <3> through a closed loop C is defined by 

and thus real valued. 

The complex conjugate of the complex number z is denoted by z. Furthermore, $t(z) 
and $s(z) denote real and the imaginary part of z, respectively. 

We consider the Dirac equation in Euclidean space-time. In D = 2 we choose the 2x2 
Dirac matrices 

71 = ( I I ) , 72 = ( ° ~' ) , 75 = -*7i72 = ( I \ y (A.67) 
In D = 4 we use the chiral representation for the Dirac matrices: 








il \ 




' a, \ 




74 = 


I -it 


o J 


, 7i = 


^< o 


, i 



1 
-1 



,(A.68) 



where 1 is the 2x2 unit matrix and <jj are the (hermitian) Pauli matrices. 

B Dirac zero modes on S 2 

We consider the Dirac operator on the sphere 5* 2 with radius R. We use stereographic 
coordinates Xu2 on S* 2 defined by 

Vi = 2 Xi, i = l,2, y 3 = p X i?, x 2 = x 2 + x 2 , (B.69) 

where y — (yi,y2, 1/3) is the vector in M 3 of the corresponding point of the embedded sphere 
S" 2 with radius R. The metric on S 2 in stereographic coordinates has the form 

ds 2 = Q 2 R {dxl + dx 2 ) , (B.70) 

where Qr — 2R 2 /(R 2 + x 2 ). The Dirac operator D in these coordinates is given by 



D = Q R 3/2 DQ][ 2 , (B.71) 

where D is the Dirac operator on R 2 . Therefore the zero modes ip of the Dirac operator 
on S 2 are related to the zero modes ip of the Dirac operator on R 2 by 

^ = n R 1/2 ^j. (B.72) 
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With z = x\ + 1x2 the zero mode of the Dirac operator in the presence of two smeared out 
center vortices at the points z = a and z = b reads 



Vi = 0, ^ 2 = - y/zz + FP_ 

j/{z - a){z -a)+ e 2 ^{z - b){z - b) + e 2 

This zero mode is normalizable with respect to the measure on S 2 which is given by 



C Uniqueness of \2 from periodicity properties 

The existence of an analytic function with the periodicity properties (|4.61|) is seen by 
choosing 

X2 (z) = 6{z + l -r- ~Q( Zl + z 2 ),ir) , (C.74) 

cf. eq. ( |4.33| ). To show the uniqueness (up to a constant factor) of X2 we assume that there 
is another analytic function Xi satisfying equation ( f4.61|) . Now consider the meromorphic 
function f(z) := X2.{ z ) / Xi{ z ) • This is an elliptic function with periods 1 and it. But f(z) 
has only a single pole in the fundamental domain (0 < Re(z) < 1 , < ^s(z) < r) at the 
zero of 9(z + | r — f^t^i + z 2 ),ir). But an elliptic function has to have at least two poles 
in the fundamental domain or it is a constant |[27|| . Hence, we infer that f(z) is a constant 
and X%( z ) is ( U P ^° a constant factor) given by ( |C.74j ). 



D Multi vortex solution on T 2 



Let us assume we have a number of thick vortices and anti-vortices at the points Zk , k = 
1, . . . , m + and z\ , I — m + + 1, . . . , m + + m~ , respectively, where the total number m + + 
vrC = m of vortices is even. The corresponding gauge potential reads 



djj2<f>& z k)- Yl <f>& z i)) > (D.75) 



z 4 



,fe=l l=m++l 



where <p(z,Zk) is defined by equation ( |4.49| ) . The fermionic zero modes ipi/^z) are given 
by 

1/4 



V>i = Yl(0t(z,z k )6+(z,z k ) +6;(z,z k )6-(z,z k ))^) x 

-1/4 

J] (e+(z, Zl )e+(z, zT) + e-(z, Zl )e-( Z , ij)) Xl (z), 

=m++l 



fe=i 

X 

i=m++l 



fe=l 
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m++m 

J] [9+(z, Zl )e+(z, Zl ) + e-(z, Zl )e-( Z , ij)) ^) 



X 

i=m++l 



where 6* is defined by equation ( J4.5U[ ) . The periodicity properties of ipx/2, cf. eq. ( 4.55| ) 



define the periodicity properties of the analytic functions Xi/2(- 2: ) 

Xl ( z + iT ) = e -((AW 2 )(-^(^/2))^(E?I^E^ ++1 .)) Xi(z)j (D . 76) 

+ = (D.77) 

X2{z + it) = ' A ^i= m ++i l JJx 2 (z), (D.78) 

X2(z + 1) = X 2(z), (D.79) 

where Am = m + — m~ and m = m + + m~. If Am > we find a Am/2-dimensional vector 
space of left-handed zero modes ip2 with analytic functions X2( z )- In the case Am < we 
find a Am/2-dimensional vector space of right-handed zero modes with analytic functions 
Xx{z). The functions X1/2 are then given by products of theta functions. In the case 
Am > we obtain 

Am/2 

X2(z) = ]Je(z + it/2 - %, it) , (D.80) 

3=1 

where the complex numbers Zj have to fulfill the conditions 

Am/2 Am/2 / m+ m \ 

»(X)f i ) = 0, 9f(X)%-)=2 9 5>- E *')■ (D.81) 

i=l i=l \fc=l Z=m++1 / 



The set of functions \2 given by equations ( P-80| , p.81|) forms an Am/2-dimensional vector 
space. 
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